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1. Introduction 

The inclusive rare B Xgj decay is a natural framework for high precision studies of FCNC, 
thanks to its low sensitivity to non-perturbative effects. As a loop induced process in the Standard 
Model (SM), it is highly sensitive to new physics In order to obtain stringent constraints on 
extensions of the SM from this decay, accurate measurments and precise theoretical predictions 
with a good control of perturbative and non-perturbative corrections have to be provided. 
On the experimental side, the latest measurements by CLEO, Belle and BaBar ^ have been com- 
bined by the Heavy Flavor Averaging Group into the current world average (WA) that reads for a 
photon energy cut of Ej> 1.6 GeV [Q] in the B-meson rest-frame 

=^(B^X,7)-;,,,^,= (3.55 ±0.24+°:?^ ±0.03) x (1.1) 

where the first error is given by the statistic and systematic uncertainty, the second one is due to the 
theory input on the shape function, and the third one is caused by the b ^ dj contamination . This 
average is in good agreement with the recent theoretical estimate including known next-to-next-to- 
leading-order (NNLO) effects [|] 

3g{B^X,y)%l,,^^^ = (3.15±0.23) x 10-^ (1.2) 

where the error consists of four types of uncertainties added in quadrature: non-perturbative (5%), 
pai^ametric (3%), higher-order (3%) and mc-interpolation ambiguity (3%). The total error of the 
present experimental WA of about 7% in Eq. ( |1.1[ ) is expected to be reduced at future B factories to 
approximately 5%. In view of this accuracy, SM calculations need to be improved with the same 
precision level by completing the NNLO QCD program. 

QCD corrections to the partonic decay rate Y{b ^ sy) contain large logarithms of the form 
a" {nib) In™ {nih/Myi/), with m < n, which should be resummed with the help of renormalization- 
group techniques. A convenient framework is an effective low-energy theory obtained from the 
SM by decoupling the heavy electroweak bosons and the top quark . The resulting effective La- 
grangian is a product of the Wilson coefficients C, (;U) with local flavor-changing operators 
up to dimension six. 

A consistent calculation of b ^ sy at the NNLO level requires three steps: /) evaluation of Q(/Xo) 
at the matching scale Hq ~ by requiring equality of Green's functions in the full and the ef- 
fective theory up to leading order in (external momenta)/Mvi7 to iff {af). All the relevant Wilson 
coefficients have already been calculated [^, ^ to this precision, by matching the four-quark op- 
erators Qi,...,Q(, and the dipole operators Q-j and at the 2- and 3-loop level respectively. //) 
calculation of the operator mixing under renormalization, by deriving the effective theory Renor- 
malization Group Equations (RGE) and evolving C,(jU) from Ho down to the low-energy scale 
Hh mh, using the anomalous-dimension matrix (ADM) to ^ [a^) . Here, the 3-loop renormaliza- 
tion in the {2i , • • • , Qe} and {Qi,Qs,} sectors was found in [||, ^, and results for the 4-loop mixing 
of 2i,...,26 into Q^ and were recently provided in [p^], thus completing the anomalous- 
dimension matrix. ///) determination of the on-shell matrix elements of the various operators at 
Hh ~ nih to [aj) . This task is not complete yet, although a number of contributions is known. 
The 2-loop matrix element of the photonic dipole operator Q-j, together with the corresponding 
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bremsstrahlung, was found in [11, O], confirmed in [13] and subsequently extended to include 



the full charm quark mass dependence in [14]. In [15], the (^a^nf) contributions were found to 
the 2-loop matrix elements of Qj and Q%, as well as to the 3 -loop matrix elements of Q\ and Q2, 
using an expansion in the quark mass ratio m^^/m^. Diagrammatically, these parts are generated by 
inserting a 1-loop quark bubble into the gluon propagator of the 2-loop Feynman diagrams. Naive 
non-abelianization (NNA) is then used to get an estimate of the complete corrections of ^ (q,^) by 
replacing Uf with — |j8o. Moreover, the contributions of the dominant operators at ^ {o^^q) to the 
photon energy spectrum have been computed in [|T6|]. 

A rather important and difficult piece that is still missing to date is the complete ^ (a^^) cal- 
culation of the matrix elements of the four-quark operators Q\ and Q2. These operators contain the 
charm quark, and the main source of uncertainty at the NLO level is related to the ambiguity asso- 



ciated to the choice of scale and scheme for [18]. As these matrix elements start contributing 



for the first time at ^ (oii ), the choice of scale and scheme for m<^. is a NNLO effect in the branching 
ratio. Therefore a calculation of {sy\Q\2\b) at iff (ocj) is crucial to reduce the overall theoretical 
uncertainty m^{B —>■ XsY) . In l|l7|], the full matrix elements of Qi and Q2 have been computed in 
the large rric limit, nic ^ m},/2. Subsequently, an interpolation in the charm quark mass has been 
done down to the physical region, under the assumption that the j3o-part is a good approximation 



at mc = . This is the source of the interpolation uncertainty mentioned below Eq. ( |1.2| ). Reduc- 
ing this uncertainty requires the evaluation of the 3-loop {sY\Q\2\b) at nic = 0, whereas removing 
it involves their calculation at the physical value of lUc, namely dealing with hundreds of 3-loop 
on-shell vertex diagrams with two scales ni}, and lUc which is a formidable task. Both of these calcu- 



lations are being pursued in [19], and we will comment on the current status in the next section. An 
important subset of diagrams contributing to the virtual 3-loop on-shell calculation of {sy\Q\2\b) 
for mc 7^ is the fermionic part which constitutes a major input both for the NNA and for the 
interpolation of the non-NNA terms between nic ^ mb/2 and ruc <mb/2, and are thus crucial for 
the accuracy of Eq. (L2). A result for these diagrams was presented in [ [l5| ] assuming that rif = 5 



massless fermions are present in the quark loop inserted into the gluon propagator. An independent 
check of this calculation as well as the validity of the massless approximation, and new results for 
the missing diagrams with heavy b and c quark loops have been recently given in jlT]]. 



2. Calculation of the matrix elements {sY\Qi^2\b) 

The iff [a}) calculation of the matrix elements {sY\Qi^2\b) is done within the framework of an 
effective theory with the Lagrangian 

b) + ^V*V,hZCi{^)Q<i^)- (2-1) 
v2 ,-=1 

Adopting the operator definitions of [120], the physical operators that are relevant for our calculation 
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Figure 1: An example graph of a cut 4-loop b-quark self-energy. 



together with the size of their Wilson coefficients read 
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Qi 



Ql,2 

Q3-6 



= {sric){cr'.b), 

e/glmb{^){sLO^''bR)F^y, 
\/g,mi,{pL){^LO^''r'bR)G\ 



'liv ) 



Claim) ~ 1 , 
|C3-6K)| <0.07 
Cjinih) ~ -0.3, 
C8K)~-0.15, 



(2.2) 



where F and V stand for various products of Dirac and color matrices. A possible way of getting 
the complete matrix elements at nic = is by interfering the operators Q[ and Q2 with the magnetic 
dipole operator Qj, then cutting the resulting 4-loop propagator diagrams in all possible ways that 
contain a photon and an s-quark in the final state. In total, 506 diagrams are generated this way each 
of which involves up to 5-particle cuts if final states with cc production are considered. A sample 
graph is shown in FIG. |l|. Since charmed hadrons in the final state are excluded experimentally, the 
^{B Xsj) does not contain contributions from reproduction. Thus, a perturbative calculation 
of b ^ X/""°"7 should be done accordingly. In order to avoid logarithmic divergences resulting 
at wic = from In terms, cuts through the c-quark loop inserted into the gluon propagator have 
been kept in our calculation. Their contribution will be subtracted at the measured value of mc 
after performing the interpolation. Moreover, since only the real part of the interference between 
the matrix elements of Q2 and Q^ contributes to the decay of — > X/""""/, we do not distinguish 
between masters that differ only in their imaginary part. This reduces the number of masters to less 
than 200. Details related to their calculation will be given elsewhere [jl9|]. 

As far as the t& [aj) 3-loop virtual correction to (sjlQialb) at me 7^ is concerned, the generated 
420 vertex diagrams have been expressed through 21231 scalar integrals that depend on the scales 
mt, and m^. With the help of Laporta algorithm [pT]], they have been subsequently reduced to 
476 masters. The latter are being evaluated using a combined approach, namely Mellin-Barnes 
technique together with differential equations solved numerically. The same techniques have been 
applied in the calculation of the fermionic contribution, therefore we refer the reader to Q for 
all related details. 

3. Results for the fermionic diagrams 

As was mentioned in the introduction, we have calculated the fermionic diagrams with three 
different quark loop insertions into the gluon propagator, namely a massless as well as heavy b and 
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Figure 2: Plots of Re(i7|g2|^)i/^'*' as function o\m^/ml with M = mf, (a) and M = mc (b) and {fih = mi,, 
Hf = 1). For comparison, we also show the M = case. 



c quark loops. Since the massless case was discussed in detail in []15|], we constrain ourself here to 
the new results related to the missing contributions from heavy loops and compare them with the 
massless approximation results. As, at ff{ajnf), the matrix elements of Qi and Q2 are related to 
each other by {sY\Qi\b) = —I / {2Nc){sY\Q2\b) , we just give results for the matrix elements of 22- 
The normalization of our amplitude is defined as follows 

{sylQilb) ffi^a^nf) = ■^mbnf{sy\Q2\b)^n'f'"usR44ub (3.1) 

where m}, denotes the b-quark pole mass, e and q are the photon polarization and momentum, 
/? = (! + 75)72 is the right handed projection operator, and nf is the number of active flavors of a 
given mass. The superscript (2) counts the powers of and M =(0, m^ormc) denotes the mass 
of the quark running in the loop inserted into the gluon propagator. The plots in FIG.§ summarize 
the outcome of our calculation. It turned out that the massless approximation overestimates the 
massive b result by a large factor, and moreover, has the opposite sign. On the other hand, less 
pronounced but non-negligible effects were observed for the massive c-quark case. 



4. Conclusions 

A complete ^ [cc^] calculation of the matrix elements (sylQi^lb) is crucial to reduce the 
overall uncertainty in the current NNLO estimate of the J^{B Xv7). This calculation is being 
pursued in |jl^. Taking new results for the complete NNLO fermionic contribution into account, 
an enhancement of 1.1% for Hh = 2.5 GeV is observed in the current estimate of the branching 
ratio. 
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